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Abstract 

We consider the Markov chain approximations for a class of singular stable-like pro- 
cesses. First we obtain properties of some Markov chains. Then we construct the 
i -^h " approximating Markov chains and give a necessary condition for the weak conver- 

gence of these chains to the singular stable- like processes. 
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1 Introduction 

o 

T— I [ 

A class of singular stable-like processes X is considered in [11]. These processes X correspond 
to the Dirichlet forms 

[ £(f,f) = f Rd f Rd (f(y)-f(x)) 2 J(x,y)m( y dy)dx, 

(L1) 

I ?={f EL\R d ): S(f \f) <oo}, 

where m(dy) is the measure on the union of coordinate axes Uf =1 lRj with Rj being the i-th 
coordinate axis of R d and m restricted to each Rj being one-dimensional Lebesgue measure 
on R. The jump kernel J(x, y) satisfies 



c{x,y) 



ifj/-xGUf =1 K<\{0}; 



J(x,y) = { \*-y\ 1+ "> 

0, otherwise, 

where c(x, y) = c(y, x) and < Ki < c(x, y) < k 2 < oo for all x and y in R d . 

In this paper, we consider the Markov chain approximations for the processes X in [11]. 
In the last few years, Markov chain approximations for symmetric Markov processes have 
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received a lot of attention. Stroock and Zheng proved the Markov chain approximations 
to symmetric diffusions in [10]. We refer to [3], [7] and [5] for Markov chain approxima- 
tions to general symmetric Markov processes and [6] for Markov chain approximations to 
non-symmetric diffusions. The Markov chain approximations to X are not considered in 
the above references since the processes X have singular jump kernels. It is natural to 
ask whether X could be approximated by Markov chains. If so, under what conditions, 
such approximation holds. The main difficulty is to get the near diagonal lower bounds in 
Proposition 3.4. We use weighted Poincare inequalities to obtain these lower bounds. This 
technique first appeared in [9], see also [8] and [12]. 

The paper is organized as follows. In Section 2 we introduce notation and define Markov 
chains related to X. In Section 3, we first construct a sequence of Markov chains. Then 
we obtain heat kernel estimates, exit time estimates and the regularity for these chains. 
In Section 4 we show the Markov chain approximations for processes X. Throughout this 
paper, if not mentioned otherwise, the letter c with or without a subscript denotes a positive 
finite constant whose exact value is unimportant and may change from line to line. 



2 Preliminaries 

Let C(-, •) : 7L d x Z d — > [0, oo) be the function satisfying 

(a) C(x,y) = C(y,x) for all x, y G Z d ; 

(b) There exist positive constants K\ and k 2 such that 

f < C(x, y) < i^fjW, i£y-xe uf =1 ZA{0}, 

\ C(x,y) = 0, otherwise, 

where Zj = Ze^ with e; being the i-th vector in M. d . 
For any x and y in Z, C(x, y) is called the conductance between x and y. Set 

G x :=^C(x,y)= C{x,x + z). 

y&1 d zeu-LjZi 

We define a symmetric Markov chain Y on Z d by 

F(Y l = y\Y o = x) = ^0-, hvx,yeZ d . 

The Markov chain Y is discrete in time and in space. We next introduce the continuous 
time version of Y. Let Y be a process that waits at a point in Z d for a length of time that is 
exponential with parameter 1, then jumps according to the jump probabilities of Y. After 
that, the process Y waits at the new point for a length of time that is exponential with 
parameter 1 and independent of what has gone before, and so on. The process Y defined 
above is the continuous time version of Y. The continuous time and continuous state process 
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closely related to both Y n and Y t is the process X corresponding to the Dirichlet form (S, J 7 ) 
in (1.1). 

Let qy (t, x, y) be the transition density of Y . Since the conductance function of Y satisfies 
conditions (Al)-(A4) in [12], Proposition 2.2 in [12] implies the following result. 

Proposition 2.1 For all x and y in 1, d , there exists a positive constant c\ such that 

q Y {t,x,y) < c 1 {r d/a A 1), for allt> 0. 



3 Heat Kernel Estimates and Regularity 

In this section, we first define a sequence of Markov chains from Y. Then we obtain heat 
kernel estimates, exit time estimates and the regularity result for these chains. 

For each p > 1, set S = p~ l Z d . For each x G S and A C S, let \i 9 x = p~ d and 

p p (A) = ^y- Define the rescaled process V as 
y eA 

V t = p-%^ for t > 0. 
We see that the Dirichlet form corresponding to V is 

f £ p (fJ) = EE (f(y) -f{x)) 2 cp{x,y), 

\ F p = {f l3(S^e): £>(/,/) <oo}, 

where C p (x,y) = p a ~ d C(px, py) for all x,y G S. 

Write p(t,x,y) for the transition density of V. Then 

P(t,x,y) = p d q Y (p a t,px,py) (3.1) 

for all x, y G Z d and t > 0. 

Proposition 3.1 For all p > 1, there exists C\ such that 

p(t,x,y) <c x t- d l a . 

Proof. This follows from (3.1) and Proposition 2.1. ■ 

For each A > 1, let V x be the process V with jumps greater than A removed. Write 
p x (t, x, y) for the transition density of the truncated process V x . The argument in the proof 
of Lemma 2.5 in [12] gives the following off-diagonal upper bound for p x (t,x,y). 

Lemma 3.2 For all t > and x,y G S, there exist c\ and c 2 such that 

P x (t,x,y) < Cl t- d/a e C2t - lx - yl/x . 
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For any set A C S, define 

T A (V) = inf {t > : V t £ A] and r A (V) = inf {t > : V t G A}. 

The upper bound in Lemma 3.2 implies the following exiting time estimates for V, whose 
proof can be found in Proposition 3.4 of [2] and Proposition 4.1 of [4]. 

Theorem 3.3 For a > and < b < 1, there exists 7 = 7(0, 6) G (0, 1) stzc/i £/ia£ for any 
R > 1 and 1 G 5, 

Recall the definition of the rescaled process V. Using Proposition 2.7, Remark 2.8 and 
Theorem 2.11 in [12], we obtain the following near diagonal lower bound for p(t,x,y). 

Proposition 3.4 There exists c > such that 

p(t,x,y)>ct- d / a , 

for all t > p~ a and \x — y\ < 2t 1//Q . 

Theorem 3.3 and the proof of Lemma 4.5 in [2] imply the following lemma. 

Lemma 3.5 Given 5 > there exists k > such that ifx, y G S, and A G S with dist(x, A) 
and dist(y,A) both larger than nt 1 / 2 , then 

W{V t = y,T A <t)< 5t~ d/a p~ d . 

Proposition 3.6 For all t > p~ a , there exist c\ > and 9 G (0, 1) such that if \x — z\, 
\y — z\ < t l / a , x,y, z G S, and r > t l l a /8, then 

W x (V t = y, T B{z r) >t)> Cl r d />- d (3.2) 

Proof. This follows easily from Proposition 3.4 and Lemma 3.5. ■ 

Remark 3.7 The above proposition still holds if we replace " \x — z\, \y — z\ < t l ^ a , x,y, z G 
S " with " \x-y\< 2t l l a , x,y G 5 " and " z " m (3.2) with " x ", respectively. 

As an application of Proposition 3.6, we have 

Corollary 3.8 For each < e < 1, there exists 9 = 9(e) G (0, 1) with the following property: 
ifx, ye S with \x-y\< t 1/a , t G [0,9 a r a ), andT C B(y,t 1/a )nS satisfies /2 p (T)t- d/a > e, 
then 

F x (V t G T and T B ( y , r ) > t) > c x e. 

In the remaining of this section we show the regularity result for V . Since V is a Hunt 
process, there is a Levy system formula for it. We refer to [4] for its proof. 
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Lemma 3.9 Let f : IR + x5x5-> IR + be a bounded measurable function vanishing on the 
diagonal. Then, for all x G S and predictable stopping time T , we have 

E * [ f(s, V s )] = E x [ f T /(*, V a , y) C P (V S , y)p d ds 

Let W t = Wq + t be a deterministic process. Then Z = (W t , V t ) is the space-time process 
on M + x S associated with V. We say that a nonnegative Borel measurable function h(t, x) 
on M + x S is parabolic in an open set B C M + x S if for all open relative compact sets 
B' C B and (t, x) G B', 

q(t,x)=E^ [h{Z T{B ,. Zs) ) 



For any to > 0, by Lemma 4.5 in [4], the function q p (t, x) = p(to — t, x, y) is parabolic in 
[0,t o )xS. 



Lemma 3.10 For each 5 G (0, 1), there exists 7 = js G (0, 1) such that for t > 0, and 
x G S, if A C Q% x, r) := [t, t + 7 r a ] x (B(x, r) n S) satisfies m ^^ r)) > 5, then 



F^ x \T A (Z)<r QPi{t>x>r) (Z))>c 1 5. 

Proof. Thanks to Corollary 3.8, this follows from using similar arguments in the proof of 
Lemma 4.7 in [2]. ■ 

Lemma 3.11 There exists a positive constant c\ such that for s > 2r and (t, x) G [0, 00) x S 

Proof. For simplicity of notation, we write r for TQ P ^,x,r)- Note that 

F& X \Z T £ Q p (t jXj s)) = ¥ X (V T i B(x,s)f]S;r < 7 r°). 

By Lemma 3.9, 

F X (V T B(x,s)nS) = E X [ f C p (y u y)p d dt] < c 2 s- a E x (r). 

\V-X\>8 

On the other hand, 

1>F X (V T £ B(x,r)C)S) =E X [ f ^ C p (V t , y) p d dt] > c 3 r- a E x (r). 

\y-x\>r 

Combining these estimates gives the required inequality. ■ 

We next derive the regularity result for V, which is also needed in the Markov chain 
approximations. 
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Theorem 3.12 There exist c > and (3 > (independent of R and p) such that for every 
bounded parabolic function q in Q p (0, Xo, 4i?), 

\q(s, x) - q(t, y)\<c UH^hRT? (\t - s\ 1/a + \y - x\f (3.3) 

holds for (s, x), (t, y) G Q p (0, x , R), where \\q\\oo,R '■= sup \q(t,y)\. In particular, 

(t,2/)6[0,7(4ii) 2 ]xS 

\p(s,x 1 ,y 1 )-p(t,x 2 ,y 2 )\ < c{tt\s)-^l a {\t-s\ l ' a + \ Xl - x 2 \ + \y1-y2\f. 

Proof. With the help of Lemmas 3.10 and 3.11, we can prove (3.3) in the same way as 
Theorem 4.9 of [2]. ■ 



4 Approximations 

In this section, we first construct the approximating Markov chains and then give a neces- 
sary condition for the weak convergence of these chains to singular stable-like processes X 
corresponding to the Dirichlet forms (£, J 7 ) in (1.1). 

For x G M. d and n G N, define 

[x]n = ([nx^/n, . . . , [nx d ]/n), S n = {[x] n :xeR d } and S' n = {[x] n : x G U^M*}. 

For any x and y in S n , let C n (x,y) be conductance on S n x S n satisfying 

j^ttz < C n (x, y) < ^^H, iiy-xe S' n - {0}; 
C n (x,y) = 0, otherwise, 

Y n the Markov chain associated with C n (x, y) and (S n , J-" n ) the Dirichlet form corresponding 
to Y n . Let p n (t, x, y) be the transition density of Y n . We can extend C n (x,y) to M. d x R d 
as follows: 

C n (x, y) = C n ([x] n , [y] n ), for x, y G R d . 

If / is a function on M. d , we define its restriction to S n by R n f(x) = f(x) for x G <S n . For 
A > 0, let be the A-resolvent for Y n and C/ A the A-resolvent for X. 



For any / and g in L 2 (S n ), set (/,#) n = Exes„ f( x )9(x) n d and ||/|| 2 , n = y/(f,f) n - 

We next prove the Markov chain approximations to singular stable-like processes X. The 
proof of the following result is similar to those in [7], [3] and references therein. 

Theorem 4.1 Suppose that for each N > 1, 

C n ([x} n , [y] n )l [N -i tN] (\x - y\) dydx J(x,y)l [N -i :N] (\x - y\) m(dy) dx 

weakly in the sense of measures as n — >■ 00. Then for each x G M d and each to > the 
F^ n -laws of {Y t n ; < t < t } converge weakly to the ¥ x -law of {X t ; < t < t } which 
corresponds to the Dirichlet form (1.1). 
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Proof. The proof will be done in several steps. 

Step 1 We show that any subsequence {rij} has a further subsequence {n Jfc } such that 
{U x . R njk f} converges uniformly on compact sets whenever / £ C*o(lR d ). 

For each x £ S n , let Q n (x) = Ylt=i[ x ij x i + V n ]- If / is a function on S n , we define 
its extension to M. d by E n f which is a Lipschitz-continuous function M. d — > M and satisfies 
conditions (a) E n f(x) = f(x) for x £ S n and (b) E n f is linear in each Q n (x). A construction 
of such function E n f is available in [2]. 

For fixed / £ C (M d ), 

||C/ n A (i^(/))||oo<||i?n(/)||oo/A<||/||oo/A. 

Therefore {U x (R n f)} is uniformly bounded, so is {E n U x R n f}. For any x and y in M. d , 
\U x R n f{[y] n )-U x R n f{[x] n )\ < f t \- xt Y^\f(z)\\p n (t,[x] n ,z)-p n (t,[y] n ,z)\n- d dt 

Jo zeS n 

POO 

+ / e~ xt J2 I f(z)\\p n (t,[x] n ,z) -p n (t,[y] n ,z)\n- d dt 



<2\\f\\ oo t + c 1 t ^ d+ ^ a \[x] n -[y] n f, 

where we used Theorem 3.12 in the last inequality and C\ is a constant independent of n. 
For any e > 0, we choose to small enough such that the first term is less than e/3. Fix such 
t , we next estimate the second term. Note that \[x] n — [y] n \ < \x — y\ + 1\fdjn. We obtain 

\[x]n-[y]nf <c 2 {\x-yf + n- p ). 

For the fixed to, there exists no £ N such that c\ t ~( d+ ^)/ a c 2 < e/3 for all n > no- 
Hence 

\U r x R n f([y] n )-U x R n f([x] n )\<e 

for all n > no and \y — x\ < 1/no- Since \[x] n — x\ < y/d/n, by the definition of E n and 
Theorem 3.12, 

\E n U x R n f(x) - U x R n f([x} n )\ < c 3 n-P, for all x £ R d . 

Therefore, for any e > 0, there exists n\ £ N such that 

\E n U x R n f(y) - E n U x R n f(x)\ < e (4.1) 

for all n > n± and \y — x\ < l/n±. This implies that {E n U x R n f} is equicontinuous on M. d . 
By the Arzela-Ascoli Theorem, any subsequence of {E n U x R n f} has a convergent further 
subsequence. Therefore, any subsequence {rij} has a further subsequence {nj k } such that 
{U x . R n . /} converges uniformly on compact sets whenever / £ Co(R d ). 

Step 2 Suppose that {n'} is a subsequence that {U x ,R n >f} converges uniformly to some 
H. We show H £ T. 

For A > 0, let u n = U x R n f. Then 

£ n (u n , u n ) = (R n f, u n ) n - X\\u n \\ 2 2 (4.2) 



Moreover, 

||Au n || 2 ,„ = ||Af/ n A i? ra /|| 2 , n < \\R n fh,n < SUp \\Rnf\\ 2 ,n < OO, 

n 

where we used lim ||-R n /||2,n = II / II 2 < 00 in the last inequality. 
Therefore the right hand side of (4.2) is bounded by 

\{R n f,u n ) n \ + X\\u\\l n < ^||i2 n /||2,«||ATx n || 2 ,„ + ^||Au n ||l >n < -supH^/lll,*- 

This implies that {£ n (u n , u n )} is uniformly bounded. 



Since u n > converges uniformly to H on B(0, N) for N > 0, by assumption, 

(H(y) — H(x)) 2 J(x, y) m(dy) dx 




N-i-<\y-x\<N 

< lim sup ^2 { u n'(y) ~ u nl (x)) 2 C nl (x,y)(n'y l - 



n x,yeS n ,,\y-x\<N 

< limsup<S n (u n ',u n >) 

n'—^oo 

< OO. 



On the other hand, 



/ H 2 (x) dx<-^ SUp \\R n f\\\ n < OO. 



>B(0,N) 

Combining these estimates and letting N — > 00, we have 

S(H,H) + \\H\\ 2 2 < 00 

and thus H G T . 



Step 3 We show lim £ n '(u n >,g) = E(H,g) for all g G 

n'—too 



Since g G Cq(R ), we can choose K large enough so that the support of g is contained in 
-8(0, A'). By the Cauchy-Schwartz inequality, 



E ( u n(v) ~ u n(x))(g(y) - g(x))C n (x, y) n 

x,yES n ,\y— x\>N 



1-d 



< (£> n ,,u n )) 1/2 ( (g(y)-g(x)) 2 c n (x, y )n- 1 - d ) 

x,yeS n ,\y—x\>N 

<2\\g\\ O0 {£ n (u n ,u n )) 1/2 ( E C n {x,y)rT 1 - d ) 1 ' 



1/2 



<c A \\g\\ooK d N- a (8 n {u n ,u n )) 
<c 5 N~ a . 



x€B(0,K)nS n \y-x\>N 
1/2 



s 



Similarly, 



E ( u n(y) ~ u n (x))(g(y) - g(x))C n (x,y)n 1_ 

x , y £ S„ , I y — x | < N ~ 1 



< 1/2 ( (9(y) ~ 9(x)) 2 C n (x, y) n~ l ~ d ) 

x,y£S n ,\y—x\<N- 1 

y - x| C n (x,y)n 

x£B(0,K)nS \y~x\<N- 1 
<C6||V^||oo^ a - 2 (^(Wn,Mn)) 1/2 

< c 7 iV a ~ 2 . 

Since if G J 7 , we can choose iV large enough such that 

(#(y) - H(x)){g(y) - g(x))J(x,y)m(dy) dx\ 

\ y -x\<tlN-i,N] 

is small. Recall that {n'} is a subsequence of {n} and U x ,R n if converges uniformly to H on 
compact sets. Therefore, 

E ( u n'(y) - u n >(x))(g(y) - g(x))C n ,(x, y) (n'y 1 - d 

x,y£S n , ,N- 1 <\y-x\<N 




~> I I ( H (y) - H(x))(g(y) - g(x))J(x,y)m(dy)dx. 

Combining these estimates gives lim £ n '(u n i,g) = S(H,g). 

n'— >oo 

Step 4 We show that E(H, g) = (/, g) - X(H, g) for all g G J= and H = U x f. 
From the above three steps, 

£(H,g) = lim E n '(u n ,,g) = lim (f,g) n > - \(u n ,,g) n > = (f,g) - X(H,g) 

for all g G Co(R d ). Note that Cg(R d ) is dense in J 7 with respect to the norm (£(•, -) + || • H 2 ,) 1 / 2 , 
see Theorem 3.9 in [11]. 

8(H,g) = (f,g)-\(H,g), for all £ G T. 

This implies that ii is the A-resolvent of / for the process corresponding to the Dirichlet 
form (SjT), that is, H = U x f. According to what we have obtained so far, we know 
that every subsequence of {U x Rnf} has a convergent further subsequence with limit U x f. 
Therefore, the whole sequence U x R n f converges to U x f whenever / G C (R d ), that is, 
lim U x R n f = U x f for / G C (R d ). 

71— >OC 

Step 5 For each / G C (R d ), we show lim F?R n f = ¥ t f. 

Using the same argument as in Step 1, we see that any sequence of {F™R n f} has a 
uniformly convergent subsequence whenever / G C (M. d ). Suppose we have a subsequence 
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{n'} such that lim F£ R n f exists. Note that 

n'— >oo 

/•OO /*oo 

U*R n f= / e~ xt ¥ n t R n fdt and t/ A / = / e" A( P ( /^. 
Jo Jo 

Using the uniqueness of Laplace transform and the fact lim U x R n f = U x f for / G Co(M. d ), 

we obtain that the whole sequence {P"i? n /} converges to F t f whenever / G C*o(M d ). 

Step 6 We show the weak convergence of the P^ n -laws of {Y t n ; < t < t } for each 
to > 0. 

It suffices to show the tightness of {Y t n ; < t < to} in the space -D[0,to] and that finite- 
dimensional distributions of {Y t n ; < t < to} converge to those of {X t ; < t < to}- Let r n 
be stopping times bounded by to and a sequence of positive real numbers converging 
to 0. Then, by Theorem 3.3 and the strong Markov property, 

F^{\YZ +Sn - YZ\ >a)= PW»(|1£ - VI > «) < V [x]n {T [[x]n ,A)(Y n ) < 7M)) < & 

for all n large enough such that 5 n < 7(0, 6). Moreover, [x] n — )■ x implies the tightness of the 
starting distributions and Theorem 3.3 implies the tightness of maxt G r 0i t ] \Y t n — Y£_\ both 
under pW«. By Theorem 1 in [1], we have the tightness of the P^ n -laws of {Y t n ; < t < t }. 
Suppose / G C (lR d ). Then, for each t G [0,t ], 

= Pt/ = lim P?i^/ = lim E ^ n i? n /(F f n ). 

n— >oo n— >oo 

Thus the one-dimensional distributions of {Y t n ; <t< to} converge to those of {X t ;0 < 
t < to}. Similarly, we can prove the finite-dimensional case using the Markov property and 
the time- homogeneity of Y n and the result in Step 5. ■ 
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